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Abstract. We study experimentally the line of a single vortex in a rotating prolate 
Bose-Einstein condensate confined in a harmonic potential. In agreement with 
predictions, we find that the vortex line is in most cases curved at the ends. We 
also present measurements of the quadrupole oscillation of the condensate in presence 
of a single vortex. A theoretical treatment to account for the short time and long time 
behaviors of these modes is developed. Finally, we use these measurements to infer the 
angular momentum per particle and relate it to the shape of the vortex line. 



1. Introduction 

The rotation of a macroscopic quantum object is a source of spectacular behaviors. 
A fascinating phenomenon is the nucleation of vortices under rotation. In superfluid 
liquid helium contained in a cylindrical bucket rotating around its axis one observes 
the nucleation of quantized vortices for a sufficiently large rotation frequency Q A 
similar phenomenon occurs in Bose-Einstein condensates confined in a rotating harmonic 
trap [2 El El El Ej- These quantized vortices in superfiuids are analogous to the quantized 
fiux lines in superconductors [?J. Among the remaining problems, is the shape of a 
vortex/fiux line and the study of its time evolution. The observation of inclined fiux 
lines in superconductors was only possible due to recent advances in electron microscopy 
[H|. In gaseous condensates, one has comparably easy access to the vortex line because 
the density of the atom cloud is low. A few disordered vortex lines have been observed by 
taking tomographic images perpendicular to the long axis of a cigar shaped condensate 
|4j. An array of more than 30 vortex lines in a pancake shaped condensate has been 
observed by transverse imaging of the whole atom cloud [9j . 

In this article we report the full length observation of a single vortex line in a cigar 
shaped condensate. We find that as a result of spontaneous symmetry breaking the line 
is generally bent [lOj. Our experimental results confirm recent predictions, in which 
the shape of the vortex line minimizing the energy of the gas was derived for a given 
rotation frequency HH DEI UHl El CHI • 

We also present measurements of the angular momentum of the rotating condensate, 
based upon surface wave spectroscopy. We develop the theoretical framework to account 
for the excitation of the quadrupole surface modes of the condensate. We then relate 
the measured angular momentum of the gas to the bending of the vortex line and its 
deviation from the center. 
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Figure 1. Image of the condensate with the probe beam propagating along the z axis, 
a) Before stirring, b) Just after stirring (r = 0). c) r = 150 ms. d) r = 500 ms. e) 
r 800 ms. f) r = 1200 ms. g) r = 1800 ms. 



2. Preparation of a single vortex 

2.1. Experimental setup 

Our ^''Rb condensate is formed by radio- frequency (rf) evaporation of 10^ magnetically 
trapped atoms in the F = tuf = 2 state. The longitudinal and transverse frequencies 
of the trap are respectively 6J^/27r = 11.8 Hz and + ujy)/4:7T = uj±/27t ^ 100 Hz (the 
X axis is vertical) . Because gravity slightly displaces the center of the trap with respect 
to the magnetic field minimum, the potential in the xy plane is not perfectly isotropic 
and we measure a 1% relative difference between Ux and Uy. 

The condensation threshold is reached at ^ 300 nK, with A^c -2x10^ atoms. We 
cool to typically T ^ 90 nK which corresponds to a condensate with A^o 5 x 10^ atoms 
and a chemical potential // ^ 70 nK. It is obtained using = uq + 10 kHz as the final 
rf, where Uq is the frequency at which the trap is emptied. During the rest of the 
experimental cycle, we maintain the evaporation rf at an adjustable level, typically 
^ = ^0 + 12 kHz, to control the temperature. 

2.2. Vortex nucleation 

Once the condensate is formed, we use an off-resonant laser beam to impose on the 
trapping potential an elliptic anisotropy in the xy plane [3j . The wavelength of the beam 
is 852 nm, its power 0.1 mW and its waist 20 /xm. Acousto-optic modulators defiect 
the position of the beam in the xy plane, thereby rotating the potential anisotropy at a 
frequency of r2/27r = 70 Hz. The potential created by the stirrer can be written : 



where the axes X, Y are deduced from the fixed x, y axes by a rotation by an angle Vtt. 
M is the atomic mass and the coefficient e measures the relative strength of the stirring 
and the magnetic potential. We use typically e ^ 4%. 

We apply this "laser stirrer" for 300 ms, during which ^ 7 vortices are nucleated 
p^. The condensate then evolves freely in the magnetic trap for an adjustable time r. 
The atom distribution at r is probed destructively by switching off the magnetic trap, 
letting the cloud expand during txoF = 25 ms and performing absorption imaging. 

The preparation of a single vortex line takes advantage of the slight static anisotropy 
of our magnetic trap, so that the angular momentum is not exactly a constant of 




(1) 
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motion. In a time r ^ 1-2 s, we observe a transition from a multi-vortex condensate 
to a condensate with a single vortex (see Fig. P). This relatively long time levels the 
fluctuations that may occur during the nucleation process. Thereby we are able to 
reproduce a condensate with a single vortex on every experimental cycle. This vortex 
line can then be studied for a time r < 10 s. 

In order to obtain an image of the full vortex line, we send two imaging beams 
aligned along the y and z directions, which probe the atom cloud simultaneously 
(Fig.[2fe^). The beams are combined onto a camera with the same magniflcation. During 
the expansion, the transverse dimensions x and y of the condensate are magnifled by 
uj±tTOF ^ 15, while the longitudinal dimension is nearly unchanged It has been 

shown theoretically that the presence of a single vortex line does not alter this expansion 
and that the coordinates of the line are scaled by the same factors ^4j. 

The transverse images in Fig. [2l3,c respectively obtained for r = 7.5 s and r = 5 s 
show the vortex line as a line of lower atomic density. Clearly this vortex line is not 
straight. It rather has the shape of a wide "U" or an unfolded "N". The fraction of N 
vs. U shaped vortices depends on temperature. Mostly U's are observed at the lowest 
temperature (T = 75 nK), while ~ half of the vortices are N-shaped for T > 100 nK. 
Fig. 1211 has been obtained with a 7- vortex lattice (one vortex at the center and the six 
other ones at the summits of an hexagon, such as in Fig.Qt). We have selected an image 
for which the lattice planes are parallel with the transverse imaging beam. The dark 
trace passing through the center of the condensate is the sum of 3 vortex lines. At the 
top and bottom, 2 lines lie on top of each other. 




Figure 2. (a) Schematic of the imaging system allowing for an observation of the 
full vortex line. The initially cigar shaped condensate is imaged after 25 ms of time 
of flight leading to the inversion of the ellipticity in the xz plane. Two beams image 
the atom cloud simultaneously along the longitudinal {z) and transverse {y) directions 
of the initial cigar, (b), (c) Transverse images of a condensate with a single vortex 
line respectively with a "U" and "N" shape. The schematic represents the form of the 
vortex line, (d) Transverse image of a condensate with 7 vortices. 
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Figure 3. (a) Schematic of the extraction of dz and d± from the vortex hne. (b) 
Evolution of dz with r. Each point corresponds to a single image, (c) Variation of dz 
with the temperature T. 



2.3. Evolution of a single vortex 

In order to give a quantitative analysis of the time evolution of a condensate with a 
single vortex, we measure the distance along the z direction between the two points 
where the vortex line leaves the condensate (Fig. Efe^). Note that the ends of the vortex 
line in the transverse image are nearly always well visible with good contrast 30%), 
even if the contrast of the line at the center of the condensate is poor (below 20%). 
Therefore this measurement can be performed on most images. Normalization by the 
length of the condensate along z leads to the quantity d^. 

In Fig. Eb, we plot dz as a function of r. Each point in the plot corresponds to a 
single image. The normalized length of the vortex line d^ decreases quasi-linearly with 
time. The graph shows that after r = 7 s the vortex is still present (in 7 s the number of 
atoms is divided by a factor 2.5). This long lifetime is reminiscent of the MIT result 0, 
where a vortex array with more than 100 vortices was produced at r = 0. The number 
of vortices was divided by 4 in ^ 5 s; however a single vortex could still be detected 
after r = 40 s. In both experiments, it is clear that the decay time of the last vortex is 
much longer than the one of the initial array. 

We have repeated this experiment for different temperatures T during the free 
evolution time r. Fig.Eb displays the measured slopes as a function of temperature. We 
observe a quasi-linear decrease of d^ for all temperatures. Our results can be compared 
to the one reported in ^Hj, where the decay rate of an array of ^ 100 vortices as a 
function of temperature was measured. In [18j, a 60% increase in temperature leads to 
an increase by a factor 17 in the decay rate. This variation is more dramatic than ours, 
indicating different decay mechanisms for a single vortex and for a large vortex array. 
In the latter case, the friction between the vortices and the thermal cloud makes the 
two components of the system stick together [19j. In the case of a single vortex, the 
most probable scenario is that the rotation of the thermal component rapidly stops due 
the residual static trap anisotropy jSIIj. The vortex line is then dragged to the edge of 
the condensate due to the friction between the two components [21j . 



Dynamics of a single vortex line in a condensate 



5 



3. Angular momentum of a condensate with a single vortex 

An important parameter for the characterization of a rotating condensate is the average 
angular momentum per particle L^. This quantity is related to the transverse quadrupole 
frequencies of the condensate uj± by the formula : 

L, = ^Mrl , (2) 

where r']_ stands for the average value of + for the condensate. The physical 
basis of ^ is clear. In absence of angular momentum along the z axis, the system is 
rotationally invariant and the two transverse quadrupole modes m = +2 and m = —2 
have the same frequencies = 6J_ (and also the same damping rates 7+ = 7_). In 
presence of a vortex the symmetry between the modes m = ±2 is broken and (|2j) relates 
the corresponding lift of degeneracy — uj^ to L^. 

The result ^ was established in [22| for a rotating condensate with an arbitrary 
number of vortices in the Thomas Fermi limit. For a condensate with a well centered 
vortex it had been obtained in [2SllSlESj- In this section we discuss the measurement 
of ijj± based on a percussional excitation of the condensate and we present the relation 
between the results of our measurements of and the shape of the vortex line. 

3.1. Measurement of uj± from the precession of the condensate axes 

We first prepare the condensate with a single vortex as described above. We then excite 
it by applying the laser potential with fixed axes {X^Y = x^y in Eq. (P)) for a short 
duration to = 0.5 ms <C ujj_^. Thereby, the wave function i/jq of the condensate is 
transformed into 

^(r) = ^o(r) + #(r) with (5^(r) = -^eMcu^x^ - y^) . (3) 

This excites a superposition of m = +2 and m = —2 modes with equal amplitudes. More 
precisely. Sip can be written in terms of the usual u± and v± functions characterizing 
the quadrupole modes m = ±2 in a Bogoliubov analysis (see e.g. j2Sl) • 

5ij{r) = (2i+(r) + ^;(r)) + (^_(r) + ^* (r)) , (4) 

where e! is a dimensionless coefl&cient proportional to the strength e and to the duration 
to of the excitation. 

We then let the cloud evolve freely in the magnetic trap for a variable time t. The 
condensate wave function becomes: 

7/;(r, t) = e-^^'/^{V^o(r) + iee'^^' (e-^^+'^+(r) + e^^+'^;(r)) 

+ ie'e-^-' (e-^^-'^_(r) + e^^-'^* (r)) } , (5) 

where is the chemical potential of the condensate. In (jSj) we added a phenomenological 
damping term for each mode. We then perform the time of fiight analysis and measure 
the size R\ and the polar angle 9 of the long axis, the size of the short axis, and 
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calculate the cloud ellipticity ( = R\/Rs- After some manipulation, we find that the 
predicted values for these quantities can be written in the following form: 

cos 26 = , ^^^^ sin 29 = , ^^^^ (6) 

and 

C - 1 « ^F\t) + G\t) , (7) 

where the functions F and G are defined by 

F{t) = e-^-^ sm{uj_t) + e"^+^ sm{uj^t) , (8) 
G{t) = e-^-^ cos(cj_t) - e-^+^ cos(u;+t) . (9) 

3.2. Various time domains for the condensate precession 

We now discuss the behavior of the precessing condensate expected from Eqs. dOHZI), 
assuming that the damping rates j± of the quadrupole modes are much smaller than 
the frequencies uj±. 

3.2.1. Short time behavior. For times short compared to 7^^, the exponential decay of 
the functions F and G can be neglected and we obtain simply: 

e = ^^^-^-)' niodulo ^ , (10) 

C - 1 oc |sin + uJ-)t/2)\ . (11) 

In this case, the lift of degeneracy 6J+ — uj^ is simply obtained from the measurement 
of the precession angle 6. The oscillation of the ellipticity is a function of the average 
frequency {ujj^ -\-uj_)/2. 

We note that combining (j2|) and ffTn|) . we obtain the result 

For very short times (i.e. uj±t <C 1) this result is actually quite general, as we show 
below. It applies to any fiuid confined in a harmonic potential with binary interactions, 
provided the fiuid is initially in a stationary state with rotational invariance, described 
by a density operator p^: 

(x^), = (2/^)0 = ri/2 (xy), = {p,Py)o = 
{xpx + Pxx)o = {ypy + Pyy)o = {xpy + yp^)o = 

Note that the system can have some angular momentum (L^ = {xpy —pyx)o ^ 0). After 
the percussional excitation by the laser potential varying as — the new density 
operator of the system is at first order in the perturbation: 



P = Po - ^2 



N 



E2 2 



(13) 
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where the sum runs over the particles of the fluid. Xi and yi denote the position 
operators of the i-th particle in the xy plane. The real coeflicient 77 is a measure of the 
strength of the percussional excitation. We now calculate the time evolution of (x^)^, 
{y'^)t and {xy)t- Using Ehrenfest theorem we obtain: 

d i ^ 1 

j,{-') = ^U[H.x^]) = -{p.x + xp^). 

1=1 

The average on the right hand side is calculated in the state and we get for short 
times : 

The calculation of {xy) is slightly more involved. It requires the value of {xpy + ypx) 
to flrst order in time and in rj. We use the fact that in equilibrium {xFy + yFx)o = 0, 
where F(r) is the force created on a given particle located in r by the — 1 remaining 
particles. We then obtain 

2r]hL,t r]hL,t'^ 
{xpy + ypx)t = , {xy)t = . (15) 

Quite remarkably the results fjl4lll5j) are independent of the binary interaction potential 
between the particles. The polar angle 9 characterizing the long axis of the condensate 
is obtained from 

which is indeed equivalent to f[T^ . The present proof, valid for any fluid conflned in a 
harmonic potential with binary interactions, holds only for t <C ujj_^. As stated above, 
(pT)|) holds also for much longer times ^ t ^ ^^^) for a pure BEC in the Thomas- 
Fermi limit [22J . 



3.2.2. Long time behavior. At longer times, the exponential decays of the amplitudes 
of the modes cannot be neglected anymore. If the two quadrupole modes decay with 
the same rate (7+ = 7_), the result ([Tn|) is valid at any time. Simply the amplitude 
( of the quadrupole oscillation becomes weaker and weaker, which makes it harder to 
detect. On the contrary, if the two rates 7+ and 7_ are different, say e.g. 7+ < 7_, one 
reaches at long time a regime where the mode m = +2 has an amplitude is much larger 
than m = —2. In this case one simply gets 9 = uj^t/2 and ( — 1 oc e~^+^ 

3.3. Experimental results 

Typical evolutions of the quantities 9 and ( are plotted in Fig. 0^- They have been 
obtained for a transverse frequency uj±/27r = 98.5 Hz and a nearly straight vortex 
(r = 2 s). We observe the two successive regimes anticipated above. For short times, one 
observes the quadrupole oscillation with precessing axes. For t > 5 ms, the precession 
rate increases and the jumps of 7r/2 in become more and more rounded. This behavior 




Figure 4. (a) Angle and ellipticity as a function of time t in presence of a 
single vortex. The value of the angular momentum per particle Lz is deduced from the 
precession of the main axes observed for small t. (b) Variation of dz and d±^ as functions 
of Lz. All measurements were binned into intervals of = 0.1 and averaged. 



indicates that the mode m = —2 decays faster than m = +2, as discussed above. The fit 
using ©-(0) gives: u;+/27r = 159.5±1.0Hz, 6J_/27r = 116.8±1.0Hz, 7+ = 19.1±2.0s-i 
and 7_ = 35.7 ± 4.0 s~^. We note that the two measured frequencies satisfy the sum 
rule uj'^ + uj'^ = 4: uj1_ [22j with a good accuracy. From the diflFerence uu^ — uj^ and the 
measured size of the condensate, we infer that the angular momentum per particle is 
Lz ^ h (to within 10%) j2Z|. The study of the difference between 7+ and 7_ is outside 
the scope of this paper and it will be discussed elsewhere [^21 • 

3.4' Relationship between L^ and the shape of the vortex line 

If we are interested only in L^, we do not have to repeat the full experimental sequence 
leading to Fig.EK- It is sufficient to perform only three experimental runs, which all start 
by preparing a condensate with a single vortex as described before. For the first two 
runs, we excite the two quadrupole surface modes m = ±2 and we probe the quadrupole 
oscillation at the two first maxima of the ellipticity ( after the fiash. From the precession 
of the condensate axes, we infer the frequency difference between the quadrupole modes, 
hence L^ [221I2Z|- For the third run, we repeat the experiment without the laser ffash 
and analyze the shape of the vortex line. As above we extract the normalized length 
dz of the vortex line along the z axis. We also measure the displacement d± of the 
axial part of the vortex line (bottom of the U) from the center of the condensate, and 
normalize it by the radius of the condensate in the xy plane (Fig. EK)- Since we have 
access only to the projection of the decentering on the xz plane, we actually measure 
d± \ cosa|, where a is the azimuthal angle of the axial part of the line. We account for 
this geometrical factor by dividing the measured displacement by (| cosa|) = 2/7r. 

Fig. HJd shows dz and d± as functions of Lz/h. For clarity we group all points into 
bins of = 0.1 h and average over dz and d±. The error bars give the statistical 
variation. The data corresponding to Lz < 0.2 h are not reproducible enough and are 
omitted. A straight {dz ^ I) and well centered {dj_ <C 1) vortex line corresponds to 
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Lz ^ Ti. The graph shows that — Lz/Ti. When decreases, we measure a decent ering 
d^ < 0.15 as long as > 0.5 h. Below = 0.5 /i, d± rises to 0.3. 

3.5. Comparison with previous theoretical works 

We now compare our experimental results with recent predictions for the shape of a 
vortex line in an inhomogeneous cigar shaped condensate [Hj Q2l E21 ^ • These 
theoretical studies consist in looking for the ground state of the condensate in a frame 
rotating at an angular frequency The general conclusion is that above a critical 
frequency Qc the ground state of the system has one or several vortices. The central 
vortex is generally bent if the trap aspect ratio uj±/uJz is large compared to 1, which 
is the case in our experiment. A precursor of this bending effect can also been found 
in [2HI l2n| in which a stability analysis of a straight vortex in an elongated condensate 
showed that some bending modes have negative energy, and are thus unstable. 

Our experimental procedure is somewhat different from the one considered in these 
theoretical studies. and Q play the role of a couple of conjugate extensive/intensive 
variables just like for instance volume and pressure. In our case no rotating anisotropy is 
imposed onto the condensate during the relevant evolution. The stirring laser has been 
switched on for a short time only at the beginning of the procedure in order to set a non- 
zero angular momentum in the system. We observe the evolution of the condensate in 
our static trap, as the angular momentum of the gas slowly decays. At first sight, neither 
Lz seems to be conserved nor Q imposed by the experimental environment. However, 
the long lifetimes observed for the vortex (up to 8 s) suggest that during an interval of a 
fraction of a second (a typical thermalization time), is almost constant. Therefore the 
description of our rotating condensate at time r should rather correspond to a system 
with a fixed than to a system rotating with a fixed The states of minimal energy 
may differ between the two descriptions jEHSOj- More precisely states with an angular 
momentum notably below h are found as saddle points of the energy functional 
calculated in the rotating frame (fixed f]), while they appear as energy minimum when 
one considers a condensate as ours, with a fixed angular momentum and where the 
Q parameter is in principle not relevant [14j. 

4. Conclusion 

We have reported the complete observation of a single quantized vortex line. The U- 
shaped vortex line that we observe is remarkably similar to those predicted and plotted 
in ^2 CSl 1^ • As pointed out in 1 the bending of the vortex line is a symmetry 
breaking effect which does not depend on the presence of a rotating anisotropy and which 
happens even in a completely symmetric setup. We have given a detailed theoretical 
account of the relation between measured condensate deformations and the quadrupole 
mode characteristics (frequency and damping rate). The quadrupole oscillations of the 
rotating condensate were used to measure the angular momentum. The shape of the 
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vortex line (bending and deviation from the center) were then related to the angular 
momentum of the system. Our results should help modelling the dissipative evolution 
of a rotating Bose-Einstein condensate. 
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